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We use gauge/string duality to model some hybrid heavy-quark potentials. The potentials under 
consideration can't be described by a single Nambu-Goto string. This is why we call them "exotic". 
For E~ , the result is in quite good agreement with lattice simulations. 
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I. 



INTRODUCTION 



Static quark potentials play an important role in the spectrum and phenomenology of both heavy mesons and 
heavy baryons. In many cases they have been computed in lattice simulations, and the results reveal a remarkable 
agreement with phenomenology. 

One of the implications of the AdS/CFT correspondence [2] is that it became a new tool to deal with strongly 
coupled gauge theories and, as a result, resumed interest in finding a string description of strong interactions. 

In this paper we continue a series of studies [3, 4] devoted to the heavy-quark potentials within a five(ten)- 
dimensional effective string theory. The model is not exactly dual to QCD, but it may be useful as a guide in 
the present period where the string dual to QCD is unavailable. In [3], the quaxk-antiquark potential was computed. 
The resulting potential is Coulomb-like at short distances and linear at long range. Subsequent work [5] has made it 
clear that the model should be taken seriously, particular in the context of consistency with the available lattice data 
as well as phenomenology. In [4], the static three-quark potential and pseudo-potential were studied. The results 
support the Y-ansatz for the baryonic area law as well as the universality of the string tension. 

The question naturally arises: What happens when the model is used for computing hybrid potentials? These poten- 
tials have been the object of many studies [1], particularly in relation to hybrid mesons and baryons [6]. Importantly, 
(exotic) hybrid mesons arc one of the goals of the GlueX experiment planned at Jefferson Lab. 

We present here an example of the hybrid potentials derived from the effective string theory in higher dimensions. 
It contains the quark-antiquark potential of [3] in an appropriate limit, but it includes as well a new object called 
" defect" . A defect is nothing but a macroscopic description of some gluonic degree of freedom in strong coupling. 

The paper is organized as follows. In section II, we discuss the fivc(ten)-dimensional effective string theory. We 
begin by summarizing the theoretical background and the results of [3] for the quark-antiquark potential. Then, we 
introduce a point-like defect to macroscopically describe a class of string excitations and, as a result, some hybrid 
potentials. We go on in Sec. Ill to discuss the lattice results as well as phenomenology. We conclude in Sec. IV with a 
brief discussion of possibilities for further study. 



First, let us set the basic framework. The background metric in question is a one-parameter deformation of the 
Euclidean AdSs space of radius TZ [3] 



II. SOME HYBRID QUARK POTENTIALS VIA GAUGE/STRING DUALITY 



A. 



General formalism 




(2.1) 



^ For a review, see [1]. 
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where dx^ = dx^ + dy^ + dz^ . 5 is a deformation parameter. In addition, we take a constant dilaton and discard other 
background fields. 

As a prelude to discussing the hybrid potentials, let us briefly consider the quark-antiquark potential (the ground 
state energy of flux tube). The potential can be determined from the expectation value of a Wilson loop. We do so 
by adopting the proposal of [7] 



{W{C))r.e-^, 



(2.2) 



where S is an area of a string worldshect bounded by a curve C at the boundary of the deformed AdS space. 

As known, if C is a rectangular loop of size R x T, then {W{C)) ~ e"^^-'^)^ as T — oo. For the background 
geometry (2.1), the potential V{R) is written in parametric form as [3] 



R{X) = 2 



dvv^e^^'-''"Hl-v^e^^^^-^"^ 



(2.3) 



V{\) = 20 



dv 

A./o ^ 



-1-v' 



+ C, 



where A is a parameter taking values in [0, 1], C is a normalization constant, and g = ^^^i 
The potential is linear at long distances 



(2.4) 



y(i?) =cri? + C + o(l), 

where a = egs is a string tension. While it is Coulomb-like at short distances 



(2.5) 



n^) = -|+C + o(l), (2.6) 



where a = {2Tr)^Q/T'^{\). 



B. Exotic hybrids 



The quark-antiquark potential corresponds to the ground state energy of a single string stretched between two 
fermionic sources. The common wisdom is that excited strings (fluxes) lead to hybrid potentials [6] . The structure of 
string excitations is quite complicated and is made up of many different kinds of elementary excitations like vibrational 
modes, loops of flux, knots, kinks, disconnected loops etc. 

Adopting the view point that some excitations remain very narrow along a string, we will model them by inserting 
local objects on a string. We will call these objects defects. The reason for this is that in the presence of defects 
string embeddings of worldsheets into spacetime are not differentiable at points where defects are located. In this 
case any use of the Nambu-Goto action alone to describe string excitations is no longer appropriate. Certainly, not 
all of excitations may be described by our construction. So we arc bound to learn something if wo succeed. 

In the case of a single defect, our construction is quite similar to that of [4] for the multi-quark potentials. So, we 
place a quark-antiquark pair at the boundary points of the five(ten)-dimensional space and consider a configuration 
in which each fermionic source is the endpoint of a fundamental string. The strings join at a defect in the interior 
as shown in Fig.l. Thus, the action has in addition to the standard Nambu-Goto actions, also a contribution arising 
from the defect. It is given by 



(2.7) 
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FIG. 1: A configuration used to calculate tiie liybrid potentials. Tiie quark and antiquark are set at a; = 7?/2 and x — —R/2, 
respectively. The defect is placed at D. 



where Si denotes the action of the string connecting the i-source with the defect. 

In the lattice formulation the hybrid potentials are extracted from generalized Wilson loops in the limit T — > C)o 
[1]. In our framework, a natural proposal for the expectation value of the generalized Wilson loop is 



(W^.o„)^e-^-», (2.8) 

where S'„i„ is an area of two string worldsheets bounded by a loop C at the boundary of AdS space and glued along 
a line in the bulk. As in [3] , in the large T limit S^i^ — TV such that a dominant contribution comes from a lateral 
surface whose area is proportional to T, while two end-surfaces only provide small corrections. 
Since we are interested in a static configuration, for the strings we take 

ti{Ti)=n, Xi{a.i) = ttiai + bi . (2.9) 
Here {Ti,ai) are worldsheet coordinates. The Nambu-Goto action of the i-string is then 



S^^TQJ^ da,w^ai + yf+rf , (2.10) 



where a prime denotes a derivative with respect to (Ji. 
The action for the defect is taken to be of the form 



5,., = rV(ro) , (2.11) 

where V can be considered as its effective potential. Unfortunately, the explicit form of V is not determined only from 
the 5-dimensional metric. It requires the knowledge of the string theory dual to QCD. We will return to this issue in 
the next sections. 

The boundary conditions on the fields are given by 



xi{Q) = -R/2, xi(l)=Xo, J/i(0)=0, yi(l)=yo, '^i(0) = 0, ri(l)=r„, (2.12) 

X2{Q) = x,, X2(l)=i?/2, y2(0)=yo, y2(l)=0, r2(0)=r„, r2(l)-0. (2.13) 

These determine the coefficients and 6; in (2.9). Thus, we have 



ai= Xo + i?/2 , hi = 



a2=^R/2-Xo, b2 = Xo 



(2.14) 
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Now, wc cxtrcmizc the total action S with respect to the worldsheet fields yiiui) and ri{ai) describing the strings 
as well as with respect to Xq, yo and Tq describing the location of the defect, with the following identifications: 
= 6y2{0) = Syo and (5ri(l) = (5r2(0) = Sro- In doing so, we use the fact that there are two symmetries which 
simplify the further analysis. 

Since the integrand in (2.10) does not depend explicitly on crj, we get the first integral of Euler-Lagrange equations 

V af + Vi + r'C 

In addition, because of translational invariance along the y-direction, there is another first integral. Combining it 
with (2.15) gives 

Pi = ^i- (2.16) 
Together with the boundary conditions these equations determine the j/j's 

2/i(o-i) = , 2/2('T2) = -yoO-2 + 2/o ■ (2.17) 
Next, we extremize the action with respect to the location of the defect. After using (2.14) and (2.17), we get 



(a^o + i?/2) /i + (xo ^ i?/2) /2 = , yo(/i+/2)=0, ri(l)/i - r^(0)/2 + 0-'V'(ro) = . (2.18) 
om this it follows that ya =0. As a result, the static configuration lies in the a;r-plane. 

Now we introduce k\ = (^ ,^0+^/2 ) ^2 = ( rI-i-xo ) " '^^^^ allows us to rewrite the first equation of (2.18) as 



^ =0. (2.19) 



Obviously, it yields the unique solution fci = ^2 = k. Combining this with (2.18), we have 



1 V 

? + (r-o) = 0. (2.20) 



If we define the first integrals /i at cti = 1 and I2 at 0-2 = and integrate over [0, 1] of dui, then by virtue of (2.19) 
we get 

'^"^»=^\Sr*»^»^""""'('-i^"''""-"')"'' p-^^' 

where A = sr^ . Obviously, it is consistent only if = 0. As a result, we end up with the most symmetric configuration. 

Now, we will compute the energy of the configuration. First, wc reduce the integrals over ai in Eq.(2.10) to that 
over Tj. This is easily done by using the first integral (2.15). Since the integral is divergent at = 0, we regularize it 
by imposing a cutoff e? Finally, the regularized expression takes the form 

Vn = V(A) + 2,^ /' 5 (l - j^.'e"'-"-') - (2.22) 



^ Importantly, in this process we use the same renormalization scheme as that for V in [3]. 
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Its e-expansion is 



Vn='^ + 0(1) . 

Subtracting the - term (quark masses) and letting e = 0, we find a finite result 



1 + k 



1-v' 



+ C, 



(2.23) 



where C is the same normalization constant as in (2.4). Note that (2.23) reduces to (2.4) at V = 0, as promised in 
the introduction.^ 

C. Concrete examples 

We will now describe a couple of concrete examples in which one can develop a level of understanding that is 
somewhat similar to that of the quark-antiquark case [3]. 

1. Model A 

First, we consider the five-dimensional geometry (2.1). In this case we specify the action S^^s as that of a point like 
particle m J ds, with m a particle mass. The latter implies that 



(2.24) 



Now (2.20) takes the form 



= =K(1-A)e- 



^/T+k 

With the form (2.25), k can be explicitly computed as a function of A. Inserted in (2.21) and (2.23), this gives 



(2.25) 



where k= 



and 



ViX) = 20 



+ C, 



(2.26) 



(2.27) 



where p{X) = 1 — (1 — A)^ e~^. A simple analysis shows that the parameter A takes values in the interval [0, 1] if 
K < 1 and [A*, 1] if k > 1. Here A* is a solution of equation p(A) = 0. 

The analysis of Eqs.(2.26) and (2.27) in the two limiting cases, long and short distances, is formally similar to that 
of [4] . From the 5d perspective, the system is prevented from getting deeper into the r-direction than l/\/s (the soft 
wall). 



3 It follows from (2.20) that fe = at V = 0. 
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At long distances, the potential V{R) is linear for any value of k 



V{R) = (tR + C + A + o{1), 



(2.28) 



with the same string tension as in (2.5). 

The point here is a finite gap between V and V 



A= ]im V{R)-V{R), with A = 2K^/Qa. 

R^oo 



(2.29) 



This is opposite to what has been found in 4d effective string theories, where there is no finite energy gap at long 
distances. For example, it is decreasing with increasing distance as A [13, 14]. Note that the gap is due to the 

effective potential (2.24) evaluated at A = 1. 

A more involved analysis shows three types of asymptotic behavior at short distances: 

For < K < 1, the potential is Coulomb-like 



V{R) 



where 



4 / 



1 3 7 
2'4'4' 



+ o{i), 



, 2F1 



1 1 3 



Note that a is positive for all allowed values of k. 
For K = 1, the potential is given by 



V{R) = C + jVR + 0{R) , 



where 7 = 2q\/V1s^. 

Finally, for larger values of k, the potential is given by 



(2.30) 



(2.31) 



(2.32) 



V{R) = V{K) + A{X,)R^ + 0{R^) . 
Here ^(A*) is given by (2.27) at A = A,. At this value of A the integral can be performed to yield 



(2.33) 



V{X^) = Y (''®'^* " + V^erfi(yA;)) + C . 
The coefficient A is given by 

^(A,)=f|s5(^V7rerf(v^)-20^e-^*) \ 
Note that A is positive for all A* of interest. 



(2.34) 



(2.35) 



2. Model B 



The above example has been formulated in a fashion that interprets the defect as a particle in 5 dimensions. 

Alternatively, one can consider a defect living in ten-dimensional space. For example, a flux loop is interpreted as a 
pair of baryon-antibaryon vertices connected by fundamental strings. According to [8], in ten-dimensional space it 
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is nothing but a branc-antibrane pair. Since the branes in question arc fivebranes, it seems natural to assume that 
•Sdot ^ T5 J (fix^/g^, with r5 a brane tension. 
For subsequent applications, we will need to know the 10-dimensional geometry. One possibility here is to take [9] 

ds^ = n^w{dt^ + dx^ + dr"^) + c-^^'-'^f^doj^duj^ . (2.36) 

This is a deformed product of AdS^ and a 5-dimensional compact space (sphere) X whose coordinates are w". The 
deformation is due to the r-dependent warp factor. 

Because the fivebranes wrap on R x X, with R along the t-axis in the deformed AdSs, we have 

V(ro) =m7^e-2^''o/^■o. (2.37) 

Note that m T5 / d^u^/g^. 

We are now ready to reproduce the potential. We begin with equation (2.20) which now takes the form 

-j^L= = .[l + AX)e-^\ (2.38) 

with K as in (2.25). As before, k can be explicitly computed as a function of A. Inserted in (2.21) and (2.23), this 
gives 



and 

V{X) = 2g 



«e-2^ 



(2.39) 



C, (2.40) 



where p{X) = 1 — (1 + 4A)^ e~^'^. Now the parameter A takes values in the interval [0, Ac] if k < 3/4e-'^/'*, [0, A^] U 
[A*, Ac] if 3/4 e^/^ < k < 1, and [A*, Ac] if k > 1. Here A* and A* are solutions of equation p(A) = and Ac is a 
solution of equation p{X) = A^e^^^"*^^. Together with (2.26) and (2.27), this is our main result. 

The above derivation has been presented in a sketchy fashion. The basic idea is to consider a defect as a loop of 
flux. Such a defect being projected onto the boundary of ten-dimensional space is tiny but it spreads well in the 
internal space X. In general, the form of V can be strongly dependent of X. In the situation considered above, we 
took ad hoc ten-dimensional geometry whose use is only siipported by phenomenological applications [9]. With this 
choice of the metric, the potential (2.37) is decreasing with r. As a result, the system is allowed to go beyond the soft 
wall in the bulk. This makes it very different from single strings [3], baryon vertices [4] and Model A, where it is not 
allowed to do so. 

Now let us analyze Model B in more detail. As in [3] , it is unclear to us how to eliminate the parameter A and find 
y as a function of r. We can, however, gain some important insights into the problem at hand from the two limiting 
cases. 

The potential of Model B also shows a linear behavior at large distances 

V{R) = aR + C + A + o{l), (2.41) 

with the same string tension as in (2.5). It may look counter intuitive, but the point is that only a finite piece of the 

infinitely long string goes beyond the wall. This piece is not able to change the leading behavior at large R. But what 
it can do is modify the size of next-to-leading order corrections. The finite gap between V and V (at large R) is now 



2V^(^^e-i-2^^ +^'^dvVt^-4e2(''^-i) - 1^ . 



(2.42) 
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Like in (2.29), the first term is a contribution of the local defect. It is nothing but a value of the effective potential 
(2.37) at A = Ac- On the other hand, the second term represents a contribution of the string piece which gets beyond 
the wall. 

The short distance structure of the solution (2.39)-(2.40) turns out to be richer than that of Model A. It also varies 

depending on a value of k. 

For < K < 3/4e^/^, the potential is Coulomb-like 

V{R) = -^ + 0{l), (2.43) 

with a given by (2.31). 

For 3/4e^/'' < k < 1, the solution (2.39)-(2.40) in fact describes three different surfaces such that the expectation 

value of the generalized Wilson loop is given by ( H^ge„ ) = Y^i=i'^i^~^^' ^ith. some weights Wi. At short distances 
the F's behave as 

'-^+0(1) ifz = l, 

Vi{R)=\v{X^) + A{\^)R^ + 0{R^) ifz = 2, (2.44) 
y{K) + A{\,)R^ + 0{R^) ifi = 3. 

Here a is given by (2.31). ^^(A*) and V(A«) are the values of the function (2.40) at A = {A*, A*}. For these values of 
A the integral can be performed to yield 



V{x) = 2Q^(Ke-^'' - + V^eA{^)) +C. 



(2.45) 



Surprisingly, the coefficients in front of R^ are expressed by the same function A as that in Model A. In the large T 
limit Vi is dominant at short distances and therefore provides V(R) = Vi{R). Note that the surfaces corresponding 
to Vi and Vz exist at short distances only. 

Finally, for larger values of k, the potential behaves as 

V{R) = y(A*) + A{K)R'^ + 0{R^) , (2.46) 
where t^(A,) and A(A*) are given by (2.45) and (2.35), respectively. 

III. PHENOMENOLOGICAL IMPLICATIONS 

Having derived the expressions for the potentials, we wish to address the issue of phenomenological implications. 
Otherwise, this will remain as an academic exercise in gauge/string duality. 

A. Comparison with lattice calculations 

The hybrid potentials are being extensively studied on the lattice [1] . They are classified by irreducible representa- 
tions of the symmetry group £>4/,. In practice, it is done by three quantum numbers. One of those is a projection of 
the total angular momentum on the axis passing through the quark and the antiquark. Our construction includes the 
system of two fundamental strings joined at the defect in the bulk such that its projection on the boundary entirely 
lies on the .x-axis. Hence it has zero angular momentum. The latter means that we should look at the so called S 
potentials. The two remaining quantum numbers correspond to symmetry (antisymmetry) under interchange of the 
string ends by a rotation by tt and by inversion in the midpoint. They are (+, — ) and {g,u), respectively. In these 
notations the usual quark-antiquark potential V is represented by S+. 
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In this section, we will compare our analytic results with the data obtained in pure SU{3) lattice gauge theory [10]. 
For V defined by (2.3) and (2.4) the fitted parameters are g, s, and C. We fit their values with the potential of 
[10]. So, we get g = 0.176, s = 0.44 GeV^, and C = 0.71 GeV.'* With these parameters fixed, the potential V has 
only a single free parameter k. We then fit it with S~ from [10]. The results are plotted in Fig. 2. We see that Model 




FIG. 2: Energy levels and E„ . The dots represent the lattice data. The curves are our predictions. Left: Model A. The 
value of K is set to 2.3. Right: Model B. The value of k is set to 2000. In both the cases we use the data and notations from 
Fig.2 of [10]. So, we have as = 0.2 fm and at ~ 0.041 fm. 

A shows a quite similar behavior to the lattice. However, for distances smaller than 0.6 fm (3 in units of Qs) the 
deviations are large that makes the model not very attractive phenomenologically. Amazingly, Model B reproduces 
the lattice data for E~ very well. 



B. More discussions on the lattice and phenomenology 

One of the predictions of our model is a finite energy gap between the long string ground state and its T,^ excited 
state. Is this reasonable? No strong observational argument seems to exclude the possibility that at distances larger 
than 1.2 fm the lattice data shown in Fig.2 are well fitted by separated parallel lines. On the other hand, one may 
reasonably argue that the range 1.2 fm < R < 2.4 fm is too small for making such a conclusion. 

Let us look at another portion of lattice data of [10]. For E~ and the energy gaps above the ground state are 
shown in Fig. 3. 

For distances smaller than 2.5 fm (12.5 in units of Os), it seems reasonable to fit the data to the linear ansatz 

^ r0.083i?K-0.89 ifn = 3, 

' [ 0.073 i?/a^ -0.88 if ?i = 4 . ^ ^ 

Here n = 3 corresponds to Yj'^ and n = 4 to E^" . Combining this with the successful parameterization of the ground 
state energy ^ 



* Interestingly, the fit from the slope of the p-meson Regge trajectories provides s = 0.45 GeV^ [11]. So, it seems that the value of s 

"slowly" depends on the flavor content. 
^ It is successful if _R exceeds 0.5 fm. So, we have to drop only a couple of points from each plot shown in Fig. 3, but that has no effect on 

our analysis. 
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FIG. 3: Energy gaps. The dots represent the lattice data. The black curves and lines are the predictions of the effective string 
models of [13] and [14]. Left: The energy gap between E~ and E^. Right: The energy gap between E~ and E^. In both the 
cases we use the data and notations from Fig.l of [10]. So, we have Us — 0.2 fm and at ~ 0.041 fm. 



12R \R 

gives 



Eo{R) = <jR + C--^+o(l] (3.2) 



E„{R) = ai? + C + A„ + ^ (Z„n - ^) + o(-) , (3.3) 

where Z3 « 0.11, Z4 « 0.12 and A3 « 0.77 GeV, A4 « 0.90 GeV. Note that the parameterization (3.3) includes that 
of the ground state with Aq = 0. 

To discuss effective string theories in a four-dimensional flat space, we must bear in mind the following: unlike 
those in a curved space (for discussion, see, e.g., [12]), it is an attempt to describe long flux tubes. So, it makes sense 
in the large R limit. 

The prediction of [13] for the string energy levels is given by 



i?„(i?)-ai? + C+^(n--) +... , (3.4) 

where the dots denote higher order terms in 1/R. For n = 0, it reduces to (3.2) that, as noted above, provides the 
successful parameterization of the ground state energy. However, it turns out to be less efficient for the excited states. 
The result is displayed in Fig. 3 by the horizontal dashed lines. Obviously, for distances smaller than 2fm the lattice 
data are very far from the prediction (3.4). 

A possible way out could be to calculate higher order corrections to (3.4). In [14], it was found 



E.a{R) = aR 



27r / 1 \ 
1 + ^^(^-12) 



C, 



(3.5) 



with (3.4) given by the first two terms in the large R expansion. But it doesn't help a lot, as seen from the black 
curves displayed in Fig. 3. 

In short, there does seem to be evidence that in the case of the S potentials the energy gaps above the ground state 
are finite for long distances, up to 2.5 fm. 

We must ask what will happen if R exceeds 2.5 fm? In [10], it is stated that for i? > 2fni the energy levels nearly 
reproduce the asymptotic -k/R gaps. We think that it is still an open question. The last point on the right of the 
plots shown in Fig. 3 may be a lattice artefact because of large errors at i? = 3fm. Obviously, in the absence of these 
points the statement is far from being indisputable. It is unfortunate that at present no new data on the S's exists. 
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Recent investigations of [15] have considered closed flux tubes of length R in SU (N) lattice gauge theories/' In this 
context, there is an interesting observation that should be mentioned in relation to the above discussion. Some states 
with = 0~ quantum numbers show large deviations from the Nambu-Goto closed string spectrum. The ansatz for 
the energy levels of these states is given by 



(3.6) 

where Eq is a ground state energy and q = {0, 1}. The value of m is fitted to m = Note that the second 

term is nothing but the energy of a free particle on a circle. 

Like in the case of the E potentials, the energy gaps in (3.6) arc finite for large R. Explicitly, AEq = m + 0(l/i?^). 
But it is not the whole story. It becomes more interesting as we compare the estimates for A3 from (3.3) to m from 
(3.6). They are given by A3 = 0.77 GeV and m = 0.78 GeV.^ This immediately raises a question whether the same 
object of gluonic nature occurred in both the cases. However, the answer might be negative because of the leading 
order corrections given by 1/R and 1/i?^ terms, respectively. 

In our discussion so far, we have not said anything about the S potentials at short distances. Here our models may 
be inappropriate because of the assumption of locality. Nevertheless, it is worth doing so to see how everything looks 
and works. 

For large k, in both models the potential after extrapolation to small R behaves as 

V{R) = Vo + AR'^ + 0(R^), (3.7) 

with Vo given by (2.34), (2.45) and A given by (2.35). This is a very satisfactory result. It is known from the effective 
string theories in flat space [14, 16], the MIT bag model [17] and lattice simulations [18].^ The constant term Vo is 
scheme dependent. It includes the common renormalization constant C, such as that in the energy of the ground 
state (2.6). This makes a quantitative comparison with others difficult. On the other hand, A is scheme independent, 
so it will be interesting to test it. 

We use the values for the parameters of Model B as those in Fig. 2 to estimate the coefficients in (3.7). Then our 
prediction for S~ is 

Vb« 1-90 GeV, ^«0.03GeV^ (3.8) 

Interestingly, the estimate for Vb is in the range of the known estimates for the lightest gluelump mass (for discussion, 
see [18]). 

IV. CONCLUDING COMMENTS 

(i) Although for the S states the results of Sec. II hold in more general circumstances, in Sec. Ill A we have restricted 
ourselves to S~ , the first excited state. The second is S~ . Our results for this case are shown in Fig.4. Just as before, 
the values of g, s, and C are fitted to the data of [10] for S+. We see that for R below about 0.6 fm the lattice data are 
well described by neither Model A nor Model B. Again, Model B performs better but still leaves room for refinement. 
A possible way to do so is to consider a couple of defects as it follows from the values of k. 

The observed agreement with the lattice data for S~ is somewhat surprising and in fact difficult to understand. 
Right now we don't have any microscopic description of the defect. So, we know nothing about its quantum numbers 
(+, — ) as well as {g, u). This poses the question of how reliable our models at short distances are. 



Eg{R)=Eo{R) 



2 , /27rgy 



° We thank M. Teper for letting us know this. 

In doing so, we use a = 0.18 GeV^. 
® For more discussion, see [1]. 
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FIG. 4: Energy levels and Eg . The dots represent the lattice data. The curves are our predictions. Left: Model A. The 
value of K is set to 2.8. Right: Model B. The value of k, is set to 3800. In both the cases we use the data and notations from 
Fig.2 of [10]. So, we have a, = 0.2 fm and at « 0.041 fm. 



(ii) There arc two important distinctions between static strings connecting the quark sources on the boundary of 
the deformed AdS space, see Fig. 5. 

First, in Models A and B the defects lead to the formation of cusps at x = 0. The cusp (deviation) angle is given 




FIG. 5: Static strings bend away from the boundary into the bulk. For models A and B, the defects are at the cusps. The 
quark and antiquark are set at x — iO.lOfm. The dashed line indicates a soft wall. 



by 




cos^ = <'v- (4.1) 



where the upper expression holds in Model A and the lower expression holds in Model B. By using (2.26) and (2.39), 
one can then write 0{R) in parametric form. In Model A the cusp angle becomes zero in the limit i? — > oo, while the 
effective potential (2.24) remains finite. 

Interestingly, there has been an extensive discussion of cusps in the context of cosmic strings [19]. In this case cusps 
can be formed when small-scale structures (kinks) propagating on the string in opposite directions collide with each 
other. 
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Second, in Model B the defect is allowed to go beyond the soft wall. This might have interesting implications for 
evaluating generalized Wilson loops at finite temperature. 

(iii) The parameterization proposed in this work can be obtained from that in [13] by a simple renormalization of 
the parameters. The transformation 



C C + A„ , n Z^n , with Aq = , (4.2) 

transforms (3.4) to (3.3). 

It is interesting to see how this can be modeled within cifFective string theories in a four-dimensional flat space. A 
good starting point is to take a string ended on two heavy quarks and then load it with point-like masses.^ On the 
other hand, it is possible to go along the way of section 2. 

For example, consider two heavy quarks separated by a distance R such that each is the endpoint of an open string. 
The strings join at a point-like defect of mass m. In the static gauge the leading terms of the derivative expansion of 
low-energy effective action are now 



S = -aR J dr-m J dr + . . . , (4.3) 

where the dots denote higher order terms. Obviously, the second term leads to a finite energy gap above the ground 
state at long distances AE = m + 0{1/R) that is equivalent to the first transformation in (4.2). The challenge is to 
see whether two-derivative terms in (4.3) lead to the second transformation. 
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